The Eckart frame Hamiltonians in the three-body problem 

A. V. Meremianirj*] 

Department of General Physics, Voronezh State University, 394006, Voronezh, Russia 

(QDated: November 19, 2012) 

Abstract 

The Eckart frame is used to separate out the collective rotations in the quantum three-body 
problem. Explicit expressions for the corresponding rotational and vibro-rotational (i.e. Coriolis) 
Hamiltonians are derived. Special attention is paid to the situation when two principal moments 
of inertia are equal in the equilibrium configuration. 
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I. INTRODUCTION 



The theoretical study of various properties of molecules and atomic nuclei requires the 
separation of their motion into the vibrational and rotational parts. Such a separation can 
be achieved by introducing the rotating reference frame (the body-frame, BF) whose axes 
are defined by the orientation of the whole system in space. As a result, the motion of 
particles with respect to BF can be considered as "vibration" while the motion of BF itself 
can be understood as the overall "rotation". At this stage the question occurs as to how 



to choose BF in an optimal way 



261 ] . An optimal BF would be that in which the couplings 



between the vibrations and the overall rotation are minimal. 

C. Eckart introduced BF 1] which minimizes the vibro-rotational (or, Coriolis) couplings 
in case when particles of the system perform small vibrations close to some equilibrium 
configuration. The Eckart BF is that in which the following relation is fulfilled 

m x [Rj x RS e<?) ] + m 2 [R 2 x R^] + m 3 [R 3 x R 3 eg) ] = 0, (1) 

where R J - e<? ' ) is the equilibrium value of the position vector Rj of the particle with the mass 



mj. The above equation is usually called "the second Eckart condition" 

a a. 

(The first 

Eckart condition is the requirement of the origin of BF to be at cm. of the system.) 

The Eckart BF has been studied for a long time 
rotational decomposition of the Hamiltonian corresponding to Eckart BF has some attractive 
features. One of which is that the Coriolis Hamiltonian in the Eckart frame is small for small- 
amplitude vibrations. Despite the broad literature on the topic some problems still remain. 
Namely, existing results related to the three-body problem are rather cumbersome. They 
also do not exhibit simple behavior in the rigid-body limit when vibrations tend to zero. In 
this limit the Coriolis (i.e. vibro-rotational) part of Eckart frame Hamiltonian vanishes and 
the rotational part of the kinetic energy becomes that of a rigid body. 

The goal of this article is to present the vibro-rotational Hamiltonians in the Eckart frame 
in most simple and compact form such that their rigid-body limit would be obvious. The 
consideration is based on the general formalism developed in 
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to 
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3M20l| . The reason is that the vibro- 



23 



24). 



The procedure of the separation of rotations and vibrations consists in the transformation 
of the total Hamiltonian into three parts 

H — H + H cor + H roU (2) 
2 



where Hq is the vibrational Hamiltonian, H rot is the rotational Hamiltonian which is 
quadratic in the components of the total angular momentum L given in BF, and H cor 
is the vibro-rotational (Coriolis) Hamiltonian which is linear both in components of L and 
the vibrational linear momenta. 



II. THE KINETIC ENERGY AND ECKART CONDITION IN TERMS OF JA- 
COBI VECTORS 



Assuming that the center of mass of the system is at the origin of the coordinate frame 
the Schrodinger equation for the system of three particles can be written as 

Id 2 Id 2 
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2/ix dv\ 2/z 2 drl, 



+ U-E) tf(ri,r 2 ) = 0. 



(3) 



where U is the potential energy operator which depends on three internal variables; ri and 
r 2 are Jacobi vectors (see fig. [1]) and the reduced masses \ii,\i 2 are defined by 



1 1 1 
— + — , — 
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Hi m 2 m 3 jj 2 mi m 2 + m 3 
The second Eckart condition in terms of Jacobi vectors reads 

Hi [pi x n] + /J 2 [p 2 x r 2 ] = 0. 



(4) 



(5) 



It is convenient to introduce the so-called "mass-scaled" Jacobi vectors by replacing 
r 1 — > ri/ y/JIi, r 2 —> r 2 / \fjI2 (and analogously for pi )2 ). 




FIG. 1: Jacobi vectors for the three-body system. CM23 is the CM of the particles m 2 and 7713. 
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FIG. 2: The geometrical meaning of the Eckart condition, a) the mutual angle between Jacobi 
vectors is 9 < vr/2, b) the mutual angle is 9 > ir/2. On each figure the areas of the dashed triangles 
are equal. 

In terms of mass-scaled Jacobi vectors the Schrodinger equation becomes 

+ A 2 2 ) + U-E) ¥(n, r 2 ) = 0, (6) 

and the second Eckart condition reads 

[pi x n] + [p 2 x r 2 ] = 0. (7) 

This condition has simple geometrical meaning. Namely, according to (J7|), the areas of two 
triangles build on pairs of vectors pi, iq and p 2 , r 2 must be equal, see Fig. [2j 

The equilibrium vectors px j2 are not moving in BF and they define the basis vectors of 
BF. Let us assume that BF basis vectors are directed along the main axes of the inertia 
tensor in the equilibrium configuration. The decomposition of the vectors px i2 over the BF 
basis vectors e 1; e 2 can be written as 

pi = X! ei + y\ e 2 , 

(8) 

p 2 = x 2 e 1 + y 2 e 2 , 

where (xk, yu) are Cartesian components of pk (n = 1, 2) in BF. Now let us introduce Eckart 
vectors f*i and f 2 by 

fi = x 1 r 1 + x 2 r 2 , 

(9) 

f2 = Vi ri + y 2 r 2 . 

Inserting ([8]) into ((?]) one arrives at the following representation of the Eckart condition 

[ei x fi] + [e 2 x f 2 ] = 0. (10) 
4 



The solution of this equation has been obtained in 23( and it reads 

(11) 



1 / , [f 2 x [f x x f 2 ] 
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1 / X ft X f 2 ]] 

e 2 — "= 12 



•F V |fi x S 

where the Eckart parameter J 7 is defined by 



2 



F = V/i+/2+2|fixf 2 |. (12) 

The orthogonality and normalization of vectors ei and e 2 as well as their compliance with 
the Eckart condition ( TlOT) can be easily proved by straightforward computations. 



The Eckart parameter J 7 can also be written as 231 ] . 



^=(e 1 -f 1 ) + (e 2 -f 2 ). (13) 

This equation clarifies the meaning of the Eckart parameter. As is seen, J 7 is the sum of 
projections of Eckart vectors on the corresponding basis vectors of Eckart BF. 

III. ECKART INTERNAL VARIABLES 

We choose the set of internal variables to be the lengths of the Eckart vectors fx , / 2 and 
cosine of their mutual angle r = cos <fi — (fi • f 2 ) / (/i/ 2 ). In order of derive the expression for 
the kinetic energy operator we note the connection between the gradient operators in terms 
of Jacobi and Eckart vectors, 

^- a =Xa k +Va h a=1 ' 2 - (14) 

From this equation one can derive the expression for the kinetic energy, 

V? + = J 2 V? + /iV^, (15) 

where Vi i2 is the gradient operator with respect to fi i2 and J 2 (ii) is the equilibrium moment 
of inertia with respect to y- (x-) axis (see Appendix iBl). 

Ii = Vi+vl, h=x\ + xl. (16) 

When deriving (fT5l) we have used the fact that the axes of BF are directed along the principal 
axes of the equilibrium inertia tensor, so that 

+ £ 2 y 2 = 0. (17) 



The connection of the Eckart internal variables with the Jacobi-bond variables follows 
from the definition (Q of vectors fi )2 and is given in Appendix |A] 

The above eq. ( II 5p allows one to derive the expression for the vibrational kinetic energy 
in terms of variables /i,/2,t. Omitting details of some routine computations which are 
analogous to those presented in |24| the final result reads 

^-±M f -wrK7! + l)^-^ + '* T) ' (18) 

where U denotes the potential energy which depends only on the internal variables. 



IV. THE CORIOLIS HAMILTON! AN 



In 



23 



24J | the general expression for the Coriolis Hamiltonian has been found which is 

3 / 3 



k=l 



(19) 



Here, denotes the projection of the total angular momentum operator on the k-th axis of 
BF (k = 1, 2, 3) and the parameters C lk and B k are determined by the differential properties 
of BF basis vectors. For the Eckart BF these parameters were calculated in 231 ] . 
In the three-body problem among parameters B k only B 3 is non-zero 23j, 



h-h (fi • fg) (T j v cot « 
B 3 = — ^— , f - - | = {h ~ h) -y 2 



(20) 



? 2 |fi x f 2 | 
where the Eckart parameter J 7 is defined by (TI2"]) . 

The parameters C 7 fc depend on the choice of the internal variables. The expression for 



C 7 fc is given in 



23|, 



c lk = - E( u % e ^ E ^ ■ v «) ^ = - E( u % E( e <? ■ ^° x Vq ]) ^ ( 2i ) 

Here, e^g is the unit totally antisymmetric tensor, U _1 is a 3 x 3-matrix whose elements 
depend on £. In the three-body problem explicit form of U _1 can be obtained using results 
of |23j, 

( |fx x f 2 | + / x 2 (f x -f 2 ) C 



U" 1 



1 



F\h x f 2 | 



(f!-f 2 ) |fiXf 2 |+/| 

Ifi x f 2 



(22) 



Noting the definition (Q of Eckart vectors, expression (fl4|) for the gradient operators V c 
and the property (fTT|) one can prove that the equation 

2 



^}p a x V a ] = h [ei x Vi] + h [e 2 x V 2 ] (23) 



a=l 



is valid. The action of the operators V on the internal variables £ gives a vector lying in the 
xy-plane of BF. This fact, together with (122]) . leads to the conclusion that among coefficients 
C y k only C 7 3 is non-zero. It is 

C 73 = -y{h{e 2 ■ Vi) - hie, ■ V 2 )) £ 7 . (24) 

This equation is the main result of the present section. 

The set of C-coefficients corresponding to Eckart internal variables £ = (/i, / 2 , cos0) 
defined in Sec. IHII can easily be obtained using (l!H|) . The result reads 

C a = -Aiei'V- r(fl ' fa) 



^ (el ' f2) - jl ^ 



C 23 = -^(ei • f 2 ) = h (25) 



C 33 = -^ (j a -j 1 + B iB^(j a ^-J 1 g) 



At this stage we can write the final expression for the Coriolis Hamiltonian (119 



H cor = -i ^§ ( cos0 f I 2 / 2 -Sr - hfi-Sr 



^ 2 V "V 9/i 9/2 



sin0— d — + (h- I 2 ) cot0). (26) 
ocosm 



+ |^/ 2 - Ji + sin 

The rigid-body limit of Coriolis Hamiltonian can be obtained using ( 1A6I) which yields 

cor w 

Thus, in the limit of zero vibrations the vibro-rotational couplings disappear. This fact is 
the inherent property of the Eckart BF which makes it so popular in molecular physics. 

The expression (126!) for the Coriolis Hamiltonian simplifies when the equilibrium inertia 
moments are equal I — I\ — L 2 , 

^-^•(-*( A &- A £) + (H) (si,i ^y < 27 » 
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V. THE ROTATIONAL HAMILTONIAN 



The rotational Hamiltonian expresses as 24j. |25] 



where, as above, = (e^ • L) denotes the z-th component o 



operator in BF and X is the effective inverse inertia tensor 



24 



(28) 



the total angular momentum 



25] . This tensor is symmetric 



T. 



X 

J 1 ' 



The general expression for the tensor X for the Eckart frame has been derived in [23]. 
However, in that paper the closed form for H rot in case of three particles was not given. 
The more detailed analysis presented in 25j makes it possible to derive a simpler form for 
the effective inertia tensor and, thereby, for the rotational Hamiltonian. In the three-body 
problem the non-zero components of the effective inertia tensor X are 25 ] 

/ 2 (cot0) 2 



In 

Xl2 









fl 


J 72 


h 


sin <j) 


A 


(" 




h 


J 72 


fx 


sin 



+ 



+ 



T 2 ' 
Ji(cot <j)) 2 

T 2 • 



cos < 



X 



(J 7 sin (j)) 2 
h + h 



^(I 1 + I 2 )sm<P + I 1 ^+I 2 j^J , 



(29) 



33 



.I 2 



The rigid-body limit of the above expressions can be easily obtained noting eqs. ( 1A6|) . This 
leads to the well-known expression for the rotational Hamiltonian for a rotating planar rigid 

system 



^1 

2/i 



^2 



T 3 
^3 



2/ 2 ' 2(/i + / 2 ) 

The equations ( 1291) simplify in practically important case when two equilibrium moments 
of inertia are equal I — I\ = I2, 

I _ / 



In 



(/a sin 0) 2: 
I cos cb 



X 



22 



-12 



X 



33 



/1/2 (sin0) 2 

Thus, the rotational Hamiltonian can be written as 



(fi sin . 

21 
~ J 2 ' 



(30) 



H, 



rot 



2 (sin . 



Li 



T 2 



cos < 



fl + fl + hh 



-L X L, 



13 
J 72 ' 



(31) 



This particularly simple expression has not yet been given in literature. 
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VI. CONCLUSION 



In the present paper the vibro-rotational decomposition of the three-body Hamiltonian 
has been analysed for the body-frame defined by the Eckart condition (J5J). The advantage 
of Eckart frame is that it minimises the vibro-rotational couplings (i.e. Coriolis part of 
the total Hamiltonian) for small amplitude vibrations. The geometrical explanation of the 
second Eckart condition (jSJ) has been given in Sec. [Til see Fig- El 

The derived expressions for the Coriolis (see eq. (126|) of Sec. IIV|) and rotational Hamil- 
tonians (eqs. f l28|) . (129]) of Sec. [V]) are written in terms of Eckart internal coordinates which 
are connected (Appendix IA1 to conventional Jacobi-bond coordinates widely used in three- 



body problem 



21 



22|. The expression for the vibrational Hamiltonian in terms of Eckart 



variables is given by eq. ffl8|) . The use of the Eckart internal variables has made it possible 
to simplify the expressions for the vibro-rotational Hamiltonians. 

It turns out that in the case of equal equilibrium inertia moments the Coriolis and ro- 
tational Hamiltonians can be written in particularly simple form, see eqs. (|27l) . ( I3TT) . and 
[9]), flAlOp . Below the sum of Coriolis and rotational Hamiltonians is written explicitly in 



terms of mass-scaled Jacobi-bond coordinates 
H, 



'■trie ~\~ H ro l L$ 



d d \ (r 2 n 

cos 9 r 2 - r iTr~ + 

or i or 2 J \r\ r 2 

Li 



(sin 9) 



d 



+ 



2 (sin i 



T 2 



dcosO 
cos 9 



rir 2 



LiL 2 + 



7-2 
^3 



(32) 



where the parameter F is defined by 

F = r\ + r\ + 2r\r 2 sin 9. 

As one can easily see, the Corioolis term as well as the non-diagonal term of the rotational 
Hamiltonian vanish in the above expression in the limit r% — > r 2 , 9 — > n/2 which is the 
equilibrium configuration with two equal principal inertia moments. At the same time, the 
rotational Hamiltonian becomes that of a rotating rigid body. 

For the sake of completeness we present here also the expression for the vibrational kinetic 
energy in terms of mass-scaled Jacobi-bond coordinates, 

2 



1/1 1 

+ 







d 



-j' 2rf dri 1 dri 2 \r\ r\) dcos9 
where U is the potential energy. 



sin i 



dcos9 



+ L 7 (ri,r 2 ,cosl 



(33) 



Appendix A: Relations between Eckart and Jacobi-bond coordinates 



Below the connections between the set of Eckart and Jacobi internal coordinates are pre- 
sented. Such connections follow from the definition of Eckart vectors flU]) and the properties 
of the principal inertia axes. Omitting details of somewhat cumbersome computations we 
present here only the final results, 

fl = j^-r (r?(/i - Pi) + r 2 2 (h - p\) ~ 2(n ■ r 2 ) ( Pl • p 2 )) . (Al) 
For the vector f 2 one has 

fl = T-V (r?(I a - Pl) + r\{I 2 ~ p\) ~ 2( ri • r 2 ) ( Pl ■ p 2 )) , (A2) 
The scalar product of Eckart vectors expresses as 

(fi ' f 2 ) = ((Pi ' P2) (ri - rl) - ( ri ■ r 2 ) (p? - pg)) . (A3) 

As is seen, the scalar product (fi • f 2 ) vanishes in the rigid-body limit when the vibration 
amplitudes tend to zero. The analysis of the rigid body limit of parameters f\ and / 2 is 
slightly more complicated. Namely, in the rigid-body limit we have r x — > pi and r 2 — > p 2 
which, noting the definitions ©, leads to 

{fi) eq = %ipi + ApI + 2x 1 x 2 (p 1 ■ p 2 ). (A4) 

Now we re- write p's in terms of their coordinates according to (jHJ). This yields 

(fi) eq = [A + + ix m + x 2 y 2 f = I\. (A5) 

Here, in the course of derivations we have utilized eqs. (fl~T|) and ([16]) . 
Thus, in the limit of zero vibrations one has 

[fl)eq = h, (f2)eq = h, 

(A6) 




Clearly, the above eqs. ( 1A1I) - ( 1A3I) are invalid when the equilibrium principal moments 
of inertia are equal. According to eq. ( 1B1I) of the next Appendix the condition 1\ = I 2 is 
met only if 

Pi = P2, e = tt/2, (A7) 
10 



where 6 e is the mutual angle in the equilibrium configuration, see Fig. [TJ Thus, the vectors 
pi and p 2 are perpendicular and, hence, we can choose them to define the Cartesian basis, 
i.e. 

e 1 = ^ e 2 = Bl. (A8) 

Pi P2 

As a consequence we have that the coordinates in eq. (jHJ) become 

Xl = Pl, V2 = P2, X2=Vl = o. 

From these identities and from (jH]) we obtain that at Ii = I 2 the Eckart vectors can be 
chosen as 

fi = pr u f 2 = pr 2 , (A9) 
where p = p\ = p 2 . The corresponding Eckart parameter (fT2|) expresses as 

jr2 = p 2 ^2 + r 2 + 2riT , 2 gin e ^ _ ^ A10 ^ 

Appendix B: The equilibrium inertia moments 

The straightforward calculation of the eigenvalues of the inertia tensor leads to the fol- 
lowing relations for the principal inertia moments of the equilibrium configuration 

h = + l -J{p\-plY + {2p 1 p 2C ose e )\ 

2 2 V ( B1 ) 

h = ~ \yJ{pl-plf + {2p lP2C ose e f. 

We remind that pi >2 are mass-scaled Jacobi vectors. To obtain expressions in terms of 
conventional Jacobi vectors one should apply the replacements pi >2 — > pi >2 ^fp x 2 to the 
above equations fIBlj) . 

As is seen from (IB II) the principal inertia moments are equal if Jacobi vectors are per- 
pendicular (i.e. when 9 e = vr/2) and their lengths satisfy the equation 

Pipl = P2p\- (B2) 



[1] C. Eckart, Phys. Rev. 47, 552 (1935). 



11 



[2] L. C. Biedenharn and J. D. Louck, Angular Momentum in Quantum Physics. Theory and 
Applications (Addison- Wesley, 1981), 



[3] 
[4] 

[5] 
[6] 

[7: 

[8 
[9 
[10 

[11 
[12 

[13 

[14 

[15 



J. K. J. Watson, Mol. Phys. 15, 479 (1968). 

A. Y. Dymarsky and K. N. Kudin, The Journal of Chemical Physics 122, 124103 (pages 2) 

(2005), URL |http : //link . aip . org/link/? JCP/122/124103/1[ 

K. L. Mardis and E. L. Sibert III, J. Chem. Phys. 106, 6618 (1997). 

G. Natanson, Molecular Physics 66, 129 (1989). 

G. A. Natanson, Chemical Physics Letters 121, 343 (1985). 

S. M. Adler-Golden and G. D. Carney, Chem. Phys. Lett. 113, 582 (1985). 

R. W. Redding and F. O. Meyer, Journal of Molecular Spectroscopy 74, 486 (1979). 

F. Jorgensen, International Journal of Quantum Chemistry 14, 55 (1978). 

F. O. Meyer and R. W. Redding, Journal of Molecular Spectroscopy 70, 410 (1978). 

O. L. Weaver, R. Y. Cusson, and L. C. Biedenharn, Ann. of Physics 102, 493 (1976). 

J. D. Louck and H. W. Galbraith, Rev. Mod. Phys. 48, 69 (1976). 

A. A. Kiselev, Optika i Spektroskopiya 24, 181 (1968), eng. trans.: Optics and Spectroscopy. 
vol.24, 1968. pp. 90-2. 

S. M. Ferigle and A. Weber, American Journal of Physics 21, 102 (1953), URL 
|http : //link . aip . org/link/?AJP/21/102/T 



[16] H. Wei, J. Chem. Phys. 118, 7202 (2003). 

[17] H. Wei, J. Chem. Phys. 118, 7208 (2003). 

[18] H. Wei and T. Carrington Jr., Chem. Phys. Lett. 287, 289 (1998). 

[19] H. Wei and T. Carrington Jr., J. Chem. Phys. 107, 9493 (1997). 

[20] H. Wei and T. Carrington Jr., J. Chem. Phys. 107, 2813 (1997). 

[21] Y. F. Smirnov and K. V. Shitikova, Sov. J. Part. Nucl. 8, 344 (1977). 

[22] R. G. Littlejohn and M. Reinsch, Rev. Mod. Phys. 69, 213 (1997). 

[23] A. V. Meremianin, J. Chem. Phys. 120, 7861 (2004). 

[24] A. V. Meremianin and J. S. Briggs, Phys. Rep. 384, 121 (2003). 

[25] A. V. Meremianin, Methods of Quantum Angular Momentum Theory in the Quantum Few- 
Body Problem (LAMBERT Academic Publishing, 2011), rus: Metody kvantovoj teorii 
uglovogo momenta v kvantovoy zadache neskolkih tel. 

[26] Since the choice of BF is not unique, its definition can be considered as the gauge convention 



12 




13 



